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Abstract 

We prove an optimal stability estimate for Electrical Impedance Tomogra- 
phy with local data, in the case when the conductivity is precisely known on 
a neighborhood of the boundary. The main novelty here is that we provide 
a rather general method which enables to obtain the Holder dependence of a 
global Dirichlet to Neumann map from a local one on a larger domain when, in 
the layer between the two boundaries, the coefficient is known. 

1 Introduction 

In this paper we consider the stability issue for Electrical Impedance Tomography 
(or, as is the same, the Calderon problem) in the case of local boundary data. As a 
consequence of the recent results of uniqueness [131 IS] > there has also been a 
number of results of stability [HI [I2j [9] , either directly associated to the conductivity 
equation setting, or associated to the Schrodinger equation setting. In general, the 
problem of the optimal rate of stability it is not yet settled. 

It is expected however that, if the unknown coefficient (either the conductivity 
or the potential) is a-priori given on a neighborhood of the boundary then the rate of 
stability is of logarithmic type with a single logarithm, which is optimal, in view of 
the well-known examples by Mandache [18] . This is in fact the result of Fathallah [9] , 
which along the lines of previous uniqueness results by Lassas, Cheney and Uhlmann 
[16j and Ammari and Uhlmann [6], proves a stability result with a single logarithm 
in the setting of the Schrodinger equation. 

Here we consider the same situation directly for the conductivity equation, see 
Theorem 12.51 in the next section, with the principal aim of providing a method of 
proof which might be adapted to other inverse boundary problems with local data 



*Work supported in part by PRIN 20089PWTPS and by the JSPS International Training Pro- 
gram (ITP). 

^Dipartimento di Matematica e Geoscienze, Universita degli Studi di Trieste, via Valerio 12/1, 
34127 Trieste, Italy. E-mail: alessang@units.it 

* Department of Mathematics, Inha University, Incheon 402-751, Korea. E-mail: kskim@inha.ac.kr 



1 



and in which the unknown parameters, or objects (such as cavities, inclusions or 
cracks) are a-priori known to be located at a positive distance from the boundary. 

The main idea in this method is that, assuming that the unknown part of a 
coefficient 7 is concentrated in a subset Dec!!, and if we fix a bigger domain D 
such that D CC D CC f2, then the full Dirichlet to Neumann map A® D associated 

to D can be determined by the local Dirichlet to Neumann map associated to a 
portion £ of d£l. In fact such a dependence has a Holder rate of stability. This is 
the essence of Theorem 13.11 below. 

In Section [5] we formulate that main assumptions and state the main Theorem 
12.51 Section [3] starts with some geometrical construction needed for the precise for- 
mulation of Theorem 13 . 1 1 which . in combination with the standard stability estimate 
with the full Dirichlet to Neumann map Theorem 13.21 enables a proof of Theorem 
12.51 The final Section [5] is devoted to the proof of Theorem 13-H this is mainly 
based on the use of singular solutions and on estimates of propagation of smallness 
for solutions of elliptic equations, some ideas are borrowed from previous work in 



Let us introduce some notation and definitions. 

Given x G R n , n > 3, we shall denote x = (x',x n ), where x' = (xi, . . . ,x n -\) G 
R n_1 , x n G R. Given x G R n , r > 0, we shall use the following notation for balls 
and cylinders. 



B r (x) = {y G R n I \y - x\ < r} , B r = B r (0) , 
B' r {x') = {y' G M"- 1 I \y' - x'\ < r} , B' r = B' r (0) , 
r a ,b(») = {y = <y, y n ) € R n I \y' - x'\ < a, \y n - x n \ < b} , T a>b = r O)6 (0) . 



We shall denote by $7 a bounded open connected subset of M n . We shall assume 
that its boundary is Lipschitz according to the following definition. 

Definition 2.1 We say that the boundary of Q is of Lipschitz class with constants 
Pq, Mq > 0, if, for any P G d£l, there exists a rigid transformation of coordinates 
under which P = and 



mum- 



2 Notation and main assumptions 



nnr^ (P) = {x = (x', Xn )eTp jLt 



x n > Z(x')} 



(2.1) 



where Z : B', 



'_po_ — >• R is a Lipschitz function satisfying 

m 



Z(0) = , 



(2.2) 
(2.3) 




2 



We shall consider an open subset £ of d£l which is accessible to measurements. 
We shall require a lower bound on its smallness as follows. 

Definition 2.2 We shall say that £ has size at least do, < do < po, if there exists 
at least one point PgS such that 

dist(P,dft\£) > d . (2.4) 

Consider a conductivity coefficient 7 defined in f2 and let us assume it is a bounded 
measurable function which satisfies the ellipticity condition 

A < 7 (x) < A -1 for all x G Q , (2.5) 

for a given positive constant A. 

Definition 2.3 We define the space of localized Dirichlet data as follows 
H l c i 2 {^) = {ip e HV 2 (dn) I supp</? C £} 

— 1/2 

and we denote by H co (£) its topological dual space. We denote with (•, •) the dual 
pairing between these two spaces, based on the standard L 2 (dQ) inner product. We 
shall denote by II • |L /rr i/2,„v ,,-1/2, the operator norm on the space of bounded 

"L(H co (2-,),H ao (2-,)) 

linear operators from Hl/, 2 (E) into H C o^ 2 (Y,). 

Now, for each ip £ i?co 2 (£), consider the weak solution u E i^ 1 (S7) to the Dirich- 
let problem 

fdiv( 7 Vu) = inn, 
\u\dn = tp on dQ . 

We introduce the local Dirichlet to Neumann map A^ as the map which asso- 

1 /2 r) 

ciates to the Dirichlet data tp € H C o (£) the boundary co- normal derivative 7^ , 

where v is the outward unit normal to d£l. More precisely, we introduce the following 
definition. 

Definition 2.4 The map 

is the operator characterized by 

(A?V, i>) = f 7 Vu ■ Vv, for every ip,ip E ^ /2 (S) , (2.7) 
Jn 

where u is the solution to the Dirichlet problem 112. 6\) and v is any function in H 1 ^) 
such that v\on = ip. 



3 



We shall consider an open subset D of f2 whose boundary is also Lipschitz with 
constants pq,Mq and which is at a given positive distance from d£l, namely we 
assume 

dist(D,dQ) > po . (2.8) 

On the unknown conductivity 7 we shall assume the following a-priori regularity 
bound 

h\\ W 2,oo {n) < E . (2.9) 

and also that it is precisely known outside D. That is, we assume that we are given a 
reference conductivity 70 which satisfies (|2.5p and (12. 9ft and the unknown 7 satisfies 

7 = 70 in ft \ I). (2.10) 



Theorem 2.5 Let ft, £ and D satisfy the above stated assumptions. If 71,72 satisfy 
(pT5|) . and (j2T0|) , £/ien we /iaue 

||7i - 72||l°o (q) < w (||A£ - A^J| £(H i/2 (s))Jf -i/ S(s)) J , 

where uo{t) is an increasing function oft>0 such that 

w(t) < C| logt|" 5 /or euery < t < e" 1 , (2.11) 

/iere C > on/y depends on the a-priori data A, i?, po; -^0) ^o> diatn(O) and on n, 
whereas 5 6 (0, 1) on/y depends on n. 

3 Proof of Theorem [23] 

Before formulating the main new tool (Theorem 13.11 below) that we shall use for the 
proof of Theorem 12.51 we need to introduce some geometrical constructions. We use 
notation and some results described in [3]. 

First we introduce an augmented domain ft by attaching to ft an open set A 
lying in its exterior and whose boundary intersects <9ft on an open portion £0 CC X. 
We refer to [31 Section 6] for details. In particular, we can choose A in such a way 
that, setting ft = ft U So U A, the following properties hold. 

There exist pi,M\ > 0, only depending on po,Mo,do, such that 

(i) ft is open, connected and has Lipschitz boundary with constants p\,M\ > 0, 

(ii) There exists Q € A such that 

B 2pi (Q) cA. 



4 



Next, if we denote, for any open set E € W 1 and h > 0, 

E h = {x G E\dist{x,dE) > h} 

we observe that there exists ho > only depending on po,Mo,do such that £1^ 
is connected for every h < ho, see for instance [H Proposition 5.5]. Note that 
B P1 (Q) Cfi ft if h< Pl . 

We introduce two domains D' , D nested as follows 

Such domains can be chosen in such a way that for suitable P2,M 2 > 0, only de- 
pending on po,Mo,do, we have 

(i) \ D' and O \ D are connected, 

(ii) D' , D have C 2 boundaries, satisfying a Lipschitz condition with constants 
P2,M 2 , 

(iii) the boundaries of D, D' , D,£l p2 have mutual distance greater than /32- 

Let us incidentally note that the set D shall be used right away in the following 
statement, Theorem 13. 11 whereas the introduction of D' shall be justified during the 
proof of the same Theorem 13.11 

Now we introduce the usual Dirichlet to Neumann map A^ D associated to the 
domain D. Namely, for a conductivity coefficient 7 satisfying (|2.5p . we consider, for 
any rj € H 1 / 2 ^!)), the solution v € H l (D) to the Dirichlet problem 

fdiv( 7 V«) =0^ in 5 , 
1 v = T] on dD , 

and we define 

Affo) = 7^ „=, (3-2) 

where f is the outward unit normal to dD. Again, A® D is identified through the 
bilinear form characterization 

(A° D V, = f l^v ■ Vw , for every n , f € # 1/2 (<9I)) , (3.3) 

where u is the solution to (|3.ip and w is any function in such that = £■ 

Theorem 3.1 Ze£ f2, £, -D and .D satisfy the above stated assumptions. 1/71,72 
satisfy (f23j) . ([2^]) and (pT0|) . i/ien we /icwe 

II A 9 -'-' _ A^ll ~ ~ ^r'llA 2 _ A E Ii' 3 <"} A\ 

H A 7i A 72 ll£(ffi/2(SD),H-i/2(SD)) ^ °H A 7i A 72 ll £(H i/2 (E)jH -i/2 (s)) ' 1 ' 

where C > 0,/3 G (0, 1) on/y depend on \, E , po, M$, do, diam(£l) and on n. 
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The proof of this theorem is the content of the next section. The other main in- 
gredient for the proof of Theorem 12.51 is the following known stability result for the 
Calderon problem with full boundary data, see pQ and also, for details, [2]. 

Theorem 3.2 Let D be as above. Suppose that 71,72 satisfy (|2.5p and (|2.9p . We 

have the following stability estimate 

Wll ~ 7S|lz,»(5) < w(||A^f - ^Wc{H^{aD),H-^{db))) ' ( 3 - 5 ) 

where u is a logarithmic modulus of continuity satisfying 

u(t) < C\ logt\~ s for every < t < e _1 , 

and C > only depends on the a-priori data A, E, po, Mq, do, diam(Q) and on n, 
whereas 5 € (0, 1) only depends on n. 

Hence, assuming Theorem 13. II proven, we can conclude as follows. 

Proof of Theorem 12.51 Let us denote ||A^ - A^ 2 1| = e. By (|3.4|) . we have 

II K db_ K dD I, _ „ <r J 

H Y Si iY 72 Wc(H 1 / 2 (dD),H- 1 / 2 (dD)) - U£ 

where, without loss of generality we may assume C > 1. If, on one hand, we have 
Ce p < e" 1 then by ([33]) 

||7i -7a||L«»(n) = hi -72|| i oo ( 5 ) < 



<C|log(Ce p )r d <C[ ) |loge 
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lo gC e \ n.__i-a 



On the other hand, if Ce^ > e 1 then, trivially, 



<A- X < 



71 ~ 72||L°°(ft) 

1 /logCeA 5 ,_i 



- x \rr) |logE 

and the thesis follows. □ 



4 Proof of Theorem 13.11 

Let us begin by observing that the reference conductivity 70 can be extended to Q, 
in such a way that the ellipticity condition ()2.5[) continues to hold in all of f2 and 
that the following Lipschitz bound holds 

IMUi.^n) <E X (4.1) 
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where E\ only depends on E, po, Mo, do. The same extension can be used for any 
7 satisfying (|2.10j) . and from now on we shall replace this assumption with the 
following one 

7 = 70 in O \ D . (4.2) 



Let 71,72 be the two conductivities appearing in the statement of Theorem 13. 1\ 
let us introduce the Green's function Gi(x,y) , i = 1, 2, for the operator for div(7jV) 
in the domain Q, that is, for any y € fi, Gi(-,y) is defined as the distributional 
solution to _ 
fdiv a; (7i(-)V a: G ! i(-,y)) = -5(- -y) in Q, , 

\Gi(;y) = onan. 

As is well-known, Gi(x, y) is symmetric, it has a singularity on the diagonal {x = y} 
of the order of \x — y\ 2 ~ n , and, away from the diagonal, it is C 1,Q -regular in each of 
the two variables x,y, moreover also the mixed derivatives V x Vj / Gj(x, y) exist and 
are locally Holder continuous away from the diagonal, see for instance [10, Theorem 
8.32]. In particular we shall make use of the following energy bound. 

Proposition 4.1 For every i/Sfi and every r > we have 

|V,.G 4 (x,y)| 2 dx<Cr 2 - n , (4.3) 



lQ\B r (y) 

where C > only depends on A and n. 

Proof. A proof can be easily obtained through a Caccioppoli type inequality and 
the well-known pointwise upper bound of the Green's function [17] , details can be 
found in p2 Proposition 3.1]. □ 
Let us fix r/i 6 C l ' a {dD), i = 1, 2 for some a 6 (0, 1) and consider vi be solutions 
to the Dirichlet problems 

fdiv( 7i V^) = in D , 
I Vi = r]i on dD . 



By Green's formula, for every x € D, we obtain 

an, 

da z 

da z 



'',(•'•')= I ~,{z)^{z)Gi{x,z) - Vi(z)ji(z)^-(x,z) 



dD 



\dvi dGi 



dD 



Note that, by the regularity assumptions on the conductivities and on dD, vi,V2 are 
differ entiable, in the classical sense, up to the boundary of D and that differentiation 
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under the integral is elementarily allowed in the above formulas. Therefore, by 
Fubini's theorem, for every x € D 

Vvi(x) ■ Vv 2 (x) = 

= / ~ „lo( z ho(w)-^-{z)—^(w)V x Gi(x,z) ■ V x G 2 (x,w)da z x a w 

JdDxdD av 

dv\ dG 2 
^lo(z)jo(w) — (z)v 2 (w)V x Gi(x,z) ■ V x - — (x,w)da z x a w 
dDxdD ov ov w 

^lo( z ho(w)v 1 (z) — (w)V x — — (x,z) -V x G 2 (x,w)da z x a w 

dDxdD 

f dG dG 

+ I _'Yo(zho(w)vi(z)v2(w)V x -^—(x,z) ■V x -^(x,w)da z x a w 

JdDxdD av z OV w 

(4.4) 

where a denotes the surface measure on dD. For any z, w £ f2 \ D let us define 

S(z,w) = [ { ll (x)- l2 {x))V x G 1 {x,z)-V x G 2 {x,w)dx . (4.5) 
Jd 

Note that, by Proposition 14.11 such integral is well defined and in fact, if z, w are 
such that dist(z,dD),dist(w,dD) > R > then we have 

\S(z,w)\ < CR 2 ~ n (4.6) 

where C > only depends on A and n. Then, by a well-known identity stemming 
from (|3.3p . we have 

{{hf - A^)vi,m) = L(m - 72)Vr;i • Vv 2 dx (4.7) 



D 



consequently, by (|43j) . ([13]) 



<(A* D - A*?)Tn,V2) = I (7i - 72)V^i • V U2 dx 



= h-I 2 -h + h (4.8) 

where we denote 

h= __jo(zho(w)^-{z)—^-(w)S(z,w)da z x a w , (4.9) 

JdDxdD av 

I 2 = _ 7o ( 2 ) 7o (^)^.( 2 ;) U2 (^) — S(z,w)dcr z x cr w , (4.10) 

JdDxdD av av w 

h= [ „Jo{z)7o(w)v 1 (z)^(w)-^-S(z,w)da z x a w , (4.11) 
JdDxdD OV z 

f d d 

h= „lo{z)lo{w)vi(z)v 2 {w) — - — S(z,w)da z x a w . (4.12) 

JdDxdD av z OV w 



For z,w G A, the Green's functions G\(-, z), G 2 (-, w) have no singularity in Q and 

1/2 

also Gi(-, z)\QQ,G 2 (-,w)\dQ G iTco (£)■ More specifically, if z,w G B pi (Q), then, by 

(3D, 

where G only depends on A, po; -&/o, do and n. Thus, recalling ([2.70 . for G 
B pi (Q), S(z,w) can be rewritten as follows 



S(z,w)= / (71 (a;) -72(a;))Va;Gfi(x,^) ■ V x G 2 (x,w)dx 
Jd 

= / (7i 0*0 -72(^))V a: G ! i(x ! 2) • V x G 2 (x,w)dx ■■ 
Jn 



Hence, 

where we denote 



<(A£-A£)Gi(.,z),G 2 (.,«0) 



|5(2,io)| < Ge for every z,w G 5 Pl (Q) (4.13) 



e - H A 7i _ ^lUw 1 / 2 ^ • (4.14) 



Si ^72ll£(H c 1 /2 (S),//- 1/2 (S)) • 

On the other hand, recalling ([4.60 . the following bound holds 

\S(z,w)\ < C for every z,w G ft \ /J 7 (4.15) 
where G > only depends on A, /9o,Mo,<io and n. Moreover, we have the following. 



Proposition 4.2 For every w G Q\D, the functions S(-, w), j^-S(-,w),i = 1, . . . , n 
are weak solutions to the elliptic equation 

div(7 Vu) = in U \ D , (4.16) 

likewise, for every z€li\D, the functions S(z, ■), -J^-S(z, -),i = 1, . . . , n are weak 
solutions to the same equation. 

Proof. It suffices to verify the weak formulation of ([4.160 with an arbitrary test 
function ip G G^°(f2 \ D). This follows in a straightforward fashion, by repeated use 
of differentiation under the integral and of Fubini's theorem. Note that use is made 
of the assumption (14.21) . □ 
Using the fact that S solves an elliptic equation in each variable, we can combine 
the smallness estimate ([4.130 with the global bound ([4.150 . so to estimate the small- 
ness of S and its derivatives on dD x dD. This task can be achieved by an estimate 
of propagation of smallness, for a general discussion on this concept we refer to [J]. 
Let us fix hi < p 2 /2 only depending on po, Mq such that (fl \ D')^ is connected. 



9 



Proposition 4.3 If v is a weak solution to (|4,16p then 

\\ v Lz((n\W) kl ) ^ C \\ v \\l\B P1 m^ v K^n\W) (4 - 17) 

where C > 0,rj G (0,1) only depend on the a-priori data A, E, po, Mq, do, diam(J7) 
and on n. 

Proof. We refer to [H Theorem 5.1], of which this proposition is just a special case. 
□ 

Applying Proposition 14.31 to v = S(-,w), for any w G B pi (Q) we obtain 

\\ s (^ w )\\m { nw )hi) < ce" for ever y ™ e B pi(Q) ( 4 - 18 ) 

and by a further application of Proposition 14.31 with respect to the w variable, we 
have 

IW.OIUfi^x^^)^^- (4-19) 

Using the elliptic equation for S(-,w) and the fact that dD is contained in (0,\D')f ll 
and at a distance greater than P2/2 from d(Q\D')h 1 , by standard interior regularity 
estimates [TOJ Theorems 8.24, 8.32] we deduce 

Wz,.)\\ 

L 2 ((U\W) h ) + W^zS{z, ■)\\ L 2((Q\-jy) h ) < Ce n for every z G 3D , (4.20) 

using now the equation for S(z, •) and its first order z-derivatives, the interior reg- 
ularity estimates give 

< Ce p (4.21) 

where we denote j3 = rj 2 . Let us now combine the above bounds with (l4.8p -( j4"TT2"1) . 
For instance, we write 

I 1 = (Af m ,g) 

where 

g(z) = (Af V2 ,S(z,-)) , 
hence, with some crude majorization, 

\h\ < C\\vi\\ H i/i(dD)\\9\\HU!i(aD) - 

<C||77i|| Hl/2(9 ^ ( ll5lli«>( 9 S) + ll v 5'll L oo( 9 5 ) J • 
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By the same reasoning and by (|4.2ip we also have 

(\\9\\ L ^(dD) + \\ V 9\\ L ^ {d D)) < C\\m\\ H i/2 {d Df P , 

and consequently 

\h \ < C , e /3 ||7?i|| Hl/2{a5) ||7? 2 || H i/2 (9j 5 ) • 
Using a similar approach for the terms in (|4.10|) - (|4.12p we arrive at 

\((Af-Af) m , m )\ < Cs^WmW^^WmW^^ , (4.22) 

for any Dirichlet data rji £ C 1,a (dD), i = 1,2, and being C 1,a (dD) dense into 
H l l 2 (dD), the proof of (|3.4p and of Theorem 13.11 is complete. 
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